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Abstract 

We analyze the possibility that bubbles of quark matter surviving the 
confinement phase transition might have become colour superconducting 
due to the enormous compression suffered by them. 

Because of the relatively high temperature of the process we compare 
the initially unpaired quark matter with the colour-flavour locked alter- 
native when the extremely large chemical potential could have increased 
the critical temperature sufficiently and find that this latter phase would 
be more stable before the bubble compression stops. 

If other physical effects have not affected completely their stability, 
these bubbles might still exist today and perhaps be observed as strangelets. 



1 Introduction 

New phases of QCD are being intensively studied with the possibility of future 
phenomenological evidence both from accelerators and astrophysical observa- 
tions. 

Apart from the ordinary hadronic phase and the expected deconfined quark- 
gluon plasma (QGP) at high temperature, interesting colour superconducting 
states of quark matter are theoretically predicted at very high pressurepP|2]. 

There is agreement that the ultimate stable phase at extremely large chem- 
ical potential /i, and consequently pressure, is the colour flavour locked one 
(CFL) where couples with zero total momentum of quarks of different colours 
and flavours u d s are paired due to the fact that mass rus is no longer relevant. 
For decreasing chemical potential, before reaching the ordinary nuclear matter, 
it is possible that intermediate phases prevail such as colour crystals 3 due to 
pairing of quarks of different momentum, or Cooper pairs of only light quarks 
u d {2SC), or even single flavour diquark condensates @|. 

*On leave of absence from Centro Atomico Bariloche, San Carlos de Bariloche, Argentina. 
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During the cooling of universe, when the temperature reached T ~ 150 MeV 
presumably a first-order phase transition occurred|Sj, even though there is de- 
bate on this point, coincident with breaking of the approximate chiral symme- 
try, passing from deconfined almost free quarks in thermodynamic equilibrium 
to quarks confined in hadrons. 

If the transition was of first order, it was produced by the expansion of 
bubbles of low-temperature phase in the sea of the high-temperature one. But 
when the coalescence velocity became smaller than that of the expansion of 
universe, a reversed situation occurred, with bubbles of the surviving quark 
matter surrounded by the more stable hadronic phasefQ. 

The negative pressure due to the energy of the false vacuum |7| inside these 
bubbles makes them contract. The chemical potential which remained constant 
in both phases during the transition begins to increase inside the bubble assum- 
ing no emission of baryons. This increase is not very relevant during the first 
stage of the compression when quarks and antiquarks annihilate with emission of 
neutrinos because of their very long mean free path. But in the final stage, when 
almost all antiquarks have disappeared due to the matter-antimatter asymmetry 
generated by earlier processes in the universe, the chemical potential increase 
makes the quark pressure rise rapidly and the bubble mass starts to increase 
too. Therefore the total pressure becomes positive and tends to decrease the 
contraction velocity. 

The interval of time of the bubble contraction turns out to be small compared 
to the age ~ 10~^ sec of the universe at the confinement transition. Therefore 
the temperature may be considered constant and must be compared with the 
critical temperature for the loss of colour superconductivity. For T ^ 150 MeV 
of all the superconducting phases only the CFL one can be probably reached 
at extremely high chemical potential. This occurs when the bubble has not yet 
stopped its contraction and from this moment the CFL free energy for variable 
number of particles can be evaluated as the sum of a free-quark contribution 
with the constraint of equal density for the 3 flavours, which is slightly larger 
than the unpaired value, minus a term due to Cooper pairs gap which makes the 
superconducting state more stable. Since this free energy y = —p, the pressure 
increases even faster and the bubble contraction is very quickly stopped. 

In Sec. 2 we give the thermodynamic expressions for massless and massive 
quarks of flavour and energy densities and of pressure as functions of temper- 
ature and chemical potential, including also a description of the CFL phase. 
Sec. 3 is devoted to the calculation of the bubble contraction through the evo- 
lution of chemical potentials which determine its mass, pressure and velocity. 
Finally in Sec. 4 we comment on physical aspects which should be also taken into 
account and on the possibility that a part of these bubbles might have survived 
till our times and could be detected by high-altitude observatories. 
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2 Thermodynamics of free and paired quarks 



For the deconfined QGP phase the high T allows the approximation of neglect- 
ing the interaction among quarks so that for each fcrmion of mass m, with 2 spin 
projections and 3 colours, the number and energy densities and the pressure are 
given by 
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in terms of a cut-off A which will not appear in the final results. 

We will consider quarks u and d with a ~ and quark s with a ~ 1 at the 
confinement transition T. 

If a; < a , which will happen for s at the beginning of bubble contraction and 
for all antiquarks because /x^ < , we may expand the denominators in powers 
of £^P£ to obtain 
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On the other hand for x > a , which is needed for u and d and occurs for s 
at the end of the bubble contraction, the integrals of Eq.(2) will be the sum of 
two parts one of the same type but with lower limit exp x to allow the expansion 
in powers of "^p'^ and another between exp a and exp x where the expansion of 
the denominator is in powers of ^^^^ . The results are 
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To analyze the relative stability of states of systems with variable number 
of particles we must compare the free energy 

J E S N 

summed over quarks and antiquarks of all flavours. 

We now turn to the possibility that quarks form Cooper pairs. 

For the case where /i >> mg and T is not too high it seems clear that di- 
quark condensates are formed due to their interaction which is attractive in the 
colour antisymmetric multiplet 3, favoured in the antisymmetric zero total spin 
and flavour antisymmetry to allow total antisymmetry of fermions, i.e. 

3 

'ZrL(p)'zfL(-p))-A5]e«^%, , (6) 

where a and i are respectively colour and flavour indices and the pair has the 
same chirality because the momenta of the partners are equal and opposite. 
Analogous expression holds for R chirality. 

In this CFL phase both the gauge SU (3)c and the approximate global chiral 
SU(3)l,b. symmetry are broken giving respectively massive gluons and mesons 
as Goldstone particles. 

The gap A can be calculated in the Nambu-Jona Lasinio approximation of 
QCD where an effective 4-fermion interaction with coupling G is considered. 

For T = the estimation|S| 

Ao - , exp (gji) (7) 

may be given with G ~ oeV'^ " 

More precise evaluationsPj are done with perturbative QCD with coupling 
9 giving 

H / 3 



where the coefficient is primarily 6 = 512 tt"*^-^^', Nf is the number of flavours 
that we take 3 and as — j:^ — — — ^ with Aqcd — 200 MeV. But additional 
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corrections to b may decrease the above expression by a factor of 5 or increase 
it by 20. 

For /X ~ 1 GeV Eq.(8) gives Aq ^ 50 MeV in agreement with the result 
with Eq.(7). Even though one cannot trust the extrapolation of the perturba- 
tive calculation down to ~ IGeV , Eq.(8) shows clearly that because of the 
asymptotic freedom limit g — > for /i — s- cx), finally Aq will decrease. 

For r > there will be a non-vanishing gap provided T < Tc — 0.57Ao. 
Therefore for very large ^ colour superconductivity might be achieved even at 
the confinement temperature T ~ IbOMeV. 

To evaluate the free energy J at T = one has to sum the expression for the 
unpaired case but with the restriction of equal Fermi momentum for all flavours, 
in order to have the same number of quarks, plus a contribution coming from 
the gap[|], i.e. 



(T=Q) 



paired 



(T=0) 



(9) 



The condition of equal Fermi momenta increases Junpair compared to the 
same expression for free quarks but the difference is compensated by the last 
term of Eq.(9) for large enough fji. 

For T > the corresponding formula will have Junpair with chemical poten- 
tials such that the densities of quarks of different flavours are equal ^01 and a 
contribution of gap which vanishes for a continuous transition|H| at T = , i.e. 
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T 
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paired unpair 

The pressure of the paired quarks will be therefore larger than that of non- 
interacting quarks for T clearly below Tc. 



3 Calculation of bubble contraction 

We will consider the force applied to the bubble, accepting that the confinement 
transition was first order, causing the increase of its relativistic momentum in 
radial direction together with release of the momentum of neutrinos emitted 
due to quark- antiquark annihilation 

F=-^f^.V^ , (11) 
dt Vn/T^ ) dt ' ^ ' 

where the rest mass M comes from quark, antiquark and false vacuum energy 
contributions. The force will result from the difference Ap between the pressure 
of false vacuum and that of quarks and antiquarks inside the bubble applied to 
its surface. 
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Taking the bubble radius r negative in order to have positive velocity v 
during contraction, from Eq.(ll) we have 

, , dM M 1 dy dP^ . , , 

Ap4..^yT^ = ,— + ^-^ + — VT^ , (12) 

where y = . 

In the first stage of the bubble compression Ap > and the rest mass 
decreases because quarks and antiquarks annihilate to mantain the density de- 
termined by thermodynamic equilibrium. On its side the momentum of the 
copiously produced neutrinos, as we will see, tends to compensate the decrease 
of the first term of RHS of Eq.(12). Therefore, the second term shows an 
increase of the velocity even though smaller than without the v term. 

On the other hand, when the bubble radius becomes small, almost all the 
antiquarks have been annihilated and the produced neutrinos momentum is 
negligible whereas the quark chemical potential raises rapidly the pressure so 
that Ap < and gives also an increase of M. As a consequence the velocity 
decreases and finally the bubble compression stops. 

In our calculation we consider the system with neutral electric charge and 
make the approximations of neglecting the presence of electrons and of surface 
effects. 

The condition of zero charge density gives the baryonic density as 



Nb Nu-Nrr I f Na - N-r N, - N- 



(13) 



V V 2\ V V 

At the beginning ji << T because with this approximation Eq.(la) gives 

^.ifr» + i?r'„ , ,14) 

SO that the contribution of quarks to baryonic density is 

-TV . (15) 



Nu-Nu _ 9.84^2 



V 

Since the entropy density is s ~ T"^, the normalized matter-antimatter asym- 
metry necessary for the primordial nucleosynthesis gives 

^ - 10-- - ^ , (16) 
s 1 

determining for the transition temperature ~ 100 MeV the consistently small 
value of chemical potential /i ~ lO^^eV. 

If the bubble suffers no evaporation it will conserve the baryonic number 
during the contraction so that, beingj^ the initial radius ~ 1 cm, according to 
Eq.(15) /i will increase reaching T for a radius ~ 10^'^ cm. 

Initially the bubble will be compressed by an "external pressure" given by 
the negative of the energy density of the false vacuum pv — ^qcd ■ when /i 
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Figure 1: Rest mass M of bubble including quark-antiquark thermodynamic 
contributions at T = 150 MeV and false vacuum energy. 

surmounts T and practically no more antiquarks survive, according to Eq.(4c) 
the pressure of the 3 flavours may be estimated by ~ jfs'A*'' which will overcome 
the external one and Ap of Eq.(12) will change its sign. 

Similarly, at the beginning when /x is very small the rest mass of the bubble 
is given by an almost constant energy density according to Eqs.(3b,4b) plus the 
constant false vacuum energy density Ag^^, so that M will decrease with V. 
But when ^ > T the term ~ /x"' in the energy density will dominate making M 
increase with the bubble contraction. 

Finally, taking the effective theory as indicative for the values of chemical 
potential achieved in the process, according to Eq.(7) when ^ ^ lOT the gap 
Ao may be sufficiently large to allow T < Tc and make the bubble colour 
superconducting before it stops. 

The numerical procedure for the evaluation of Eq.(12) has been as follows. 
We started with a radius variable tq = —1 cm and zero velocity vq = 0. For the 
different steps the conservation of baryonic number Nb gave through Eq.(13) 
the chemical potential /i„ and a relation between /i^ and /ig. These two latter 
chemical potentials were determined by minimizing the energy at each radius. 
Then the pressure given by quarks and antiquarks could be calculated. The 
precise value of the negative pressure of false vacuum was adjusted to be slightly 
larger than that of quarks and antiquarks in the initial stage in order to have 
the minimal compression condition. 

Regarding the inclusion of the momentum of neutrinos in Eq.(ll) it is con- 
venient to consider the two stages of /U << T and /U ~ T. 

In the former we may approximate the energy density as a constant p = 
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Figure 2: Bubble chemical potential and internal pressure. Continuous line 
logHulMeV], broken line log(^). 



Pv + Pq,q where pq^q is due to all quarks and antiquarks, so that ^ = —p in r^. 
In the rest frame Pq^g is transformed into neutrinos which will be isotropically 
distributed with an average momentum k such that 

dEv= Pq^qA-K r'^vdt = A'Kk^a , (17) 

where (7 is surface energy density in space k. 

Taking the motion in the direction 3, the momentum seen from lab frame 
will be ^3 = ^i^^2^*^ • Therefore there will be a total momentum in lab system 

dPa = k^ [ (y) k^ sin ed9dif= , ^ ^ Pq,q A-nr'^dt , (18) 
J ^kJ vl — 

and the neutrino momentum will cancel the contribution of Pq^q in ^ of 
Eq.(12), leaving 

/^p^l-y = -pvy+ \ \r\ P-r^T- ' ^l^) 

6 1 — y dr 

In the small radius stage when x ~ 1, the antiquarks have been almost 
all annihilated so that Pg « Pq and only 2pq will transform into neutrinos. 
Therefore Eq.(12) must be used with 

dP 

^/^^^=47rr2y2p^ . (20) 

In this last part of bubble compression, when T < Tc the chemical potentials 
Pd and Ps will be determined by the condition of having equal density of u, d and 
s. The pressure of quark matter will be increased by the gap term of Eq.(lO). 
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Figure 3: Velocity of contraction of surviving quark matter bubble at T = 150 Me V. 

In this way the properties of the bubble were calculated, i.e. the rest mass 
M (Fig.l) the chemical potential and the internal pressure p (Fig. 2) and the 
velocity of contraction v (Fig. 3). 

From them it is seen that the bubble stops its compression for a final radius 
Rf slightly below 10~^cto, its rest mass starts to increase for radius slightly 
below 10~'^cm when also the increase of the internal pressure overcomes the 
external one, and in between these sizes reaches a value ~ l.bGeV so that 
T < Tc allowing the possibility that the bubble quark matter becomes colour 
superconducting. 

Finally, the interval of time for the bubble contraction 



can be calculated numerically giving At ^ 10 sec much smaller than 10 ^ sec 
so that the approximation of T ~ constant is consistent. 

4 Conclusions 

According to our model, the bubbles of quark matter surviving the confinement 
transition may reach such high values of chemical potential, larger than those in 
neutron stars, that could enter the colour superconducting CFL phase before 
their contraction stops becoming therefore stable. A more accurate analysis of 
the order of the transition between QGP and CFL phases JT] should be done 
to establish with certainty the dependence of Tc on fi. 




(21) 
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There are several physical considerations that may affect the above result. 
One is the emission of neutrons during the bubble compression which would 
make the increase of chemical potential softer. Another thing is the possible 
inclusion of a density of electrons which would alter the quark chemical poten- 
tials to satisfy charge neutrality. Finally, taking into account surface effects the 
condition of charge neutrality of the bubble may not be realistic. 

In case that the bubbles, perhaps smaller, can nevertheless become colour 
superconducting at the relatively large temperature ^ 150AIeV, one has to 
take into account that if their energy per baryon is larger than 940 MeV they 
will emit neutrons during the cooling till the present time. If this process leaves 
bubbles of atomic number A > 10^ at T ~ they would be absolutely stablefT^ 
and might be detected as strangelets by high-altitude observatories since the 
collisions in the upper part of the atmosphere would not reduce substantially 
their mass maintaining therefore their stability. 
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